The Force- Velocity Relation for Growing Biopolymers 
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The process of force generation by the growth of biopolymers is simulated via a Langevin-dynamics approach. The interaction 
forces are taken to have simple forms that favor the growth of straight fibers from solution. The force-velocity relation is 
obtained from the simulations for two versions of the monomer-monomer force field. It is found that the growth rate drops 
off more rapidly with applied force than expected from the simplest theories based on thermal motion of the obstacle. The 
discrepancies amount to a factor of three or more when the applied force exceeds 2.5 kT/a, where a is the step size for the 
polymer growth. These results are explained on the basis of restricted diffusion of monomers near the fiber tip. It is also found 
that the mobility of the obstacle has little effect on the growth rate, over a broad range. 
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PACS numbers: 87.15. Rn, 87.16. Ac, 87.17.Jj 
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I. INTRODUCTION 

The growth of biopolymers is a key ingredient in the 
crawling motion and internal transport processes of al- 
most all eukaryotic cells. They crawl among each other 
and over substrates by motion of the cytoplasm into pro- 
trusions known as lamcllipodia, filopodia, or microspikes 
according to their shapes. The forces driving the exten- 
sion of these profusions are believed to comes from the 
growth of a collection of fibers assembled from monomers 
of the protein actin. The actin fibers are approximately 
7 nm in diameter. With no opposing force, they can grow 
at velocities [El of over 1 /im/sec at physiological actin 
concentrations [^]||] of 10-50/iM; the velocities of the cell 
protrusions are typically in the range of 0.1 /im/sec. 
Actin fiber growth also can power the motion of bacte- 
ria and viruses through the cell cytoplasm. The veloci- 
ties usually range from 0.02 to 0.2 /im/sec, but velocities 
up to 1.5 /im/sec have been observed. As they move, 
they leave behind "comet tails" made up of actin fibers 
|pL»). Recent experiments have studied the minimal in- 
gredients necessary for such propulsion. For example, 
Ref. H shows that polystyrene beads coated with a cat- 
alytic agent for actin polymerization spontaneously move 
in cell extracts at velocities of 0.01 to 0.05 /im/sec, form- 
ing comet tails similar to those caused by bacteria and 
viruses. It has also been shown recently that Listeria 
and Shigella bacteria can move in a medium much sim- 
pler than a cell extract, containing in addition to actin 
monomers only the proteins Arp2/3 complex, actin de- 
polymerizing factor, and capping protein. In particular, 
myosin-type motors are not necessary for motion driven 
by actin polymerization. The minimal ingredients lead 
to velocities of 0.01 to 0.02 /im/sec; supplementation 
of this mix with other ingredients including profilin, a- 
actinin, and the VASP protein increases the velocities 



up to 0.05 /im/sec. To our knowledge, there have been 
no measurements of the force-velocity relation for grow- 
ing actin filaments. However, recent measurements of 
the actin fiber density |j] and Young's modulus (!(]] at 
the leading edge of lamellipodia would suggest forces on 
the order of 1 pN per fiber if all fibers are contributing 
equally; this is roughly equal to the basic force unit for 
fiber growth, kT/a, where k is Boltzmann's constant, T 
is temperature, and a — 2.7 nm is the incremental fiber 
length per added monomer. 

Microtubules, which are thicker fibers (22 nm) assem- 
bled from tubulin subunits, also exert forces when they 
grow. Microtubule assembly and disassembly is crucial 
in intracellular processes such as mitosis, the formation 
of cilia and flagella, and the transport of nutrients across 
the cell. Recent measurements [111 on microtubules in 
vitro have yielded explicit force-velocity curves. At zero 
force, the velocity is about 0.02 /im/sec; with increasing 
force, the velocity drops off roughly exponentially. 

It is clear that growth of the fiber against a force results 
in a lowering of the system's free energy if the opposing 
force is sufficiently small, since the exothermic contribu- 
tion from the attachment of monomers at the end of the 
polymer will outweigh the work done to move the ob- 
stacle against the external force. The critical force at 
which polymerization stops is determined by the balance 
of these two contributions. However, it is not yet under- 
stood in detail what factors determine the rate of growth 
and the maximum force at which a useful speed can be 
obtained. The basic difficulty of the polymer's growth 
process is that when the obstacle impinges directly on the 
fiber tip, there is not enough room for a new monomer 
to move in. Thus the rate of growth must be connected 
to the fluctuations of either the obstacle or the filament 
tip, which create temporary gaps between the tip and 
the obstacle. This effect has been treated explicitly in 
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the "thermal ratchet" model [Q . In this model, one as- 
sumes that the obstacle must be a critical distance a from 
the tip for growth to occur. The fiber is assumed to be 
rigid. The growth rate is obtained by solution of a drift- 
diffusion type equation. For conditions of slow growth, 
i.e. in which the time to add a monomer is much longer 
than the time it takes the obstacle to diffuse a distance 
a, this equation can be solved analytically. The forward 
growth rate is proportional to the probability that the 
obstacle-tip separation exceeds a. If depolymerization is 
sufficiently slow to be ignored, this yields the following 
dependence of the velocity v on the opposing force F: 



v oc exp (—Fa/kT) 



(1) 



where k is Boltzmann's constant and T is the tempera- 
ture. This result is equivalent to application of the prin- 
ciple of detailed balance jl3), on the assumption that 
the depolymerization rate is independent of the oppos- 
ing force. This work has been extended to flexible fibers 
at non-perpendicular incidence [p^ , p^[ , and to interacting 
systems of fibers Jlfj]. For flexible fibers, it is again found 
that the velocity is proportional to the probability form- 
ing of a gap large enough to admit the next monomer. 

It is the purpose of this paper to evaluate the force- 
velocity relation for growing fibers using a model more 
realistic than those used previously. The model used to 
derive Eq. (|l|) does not explicitly treat the diffusion of 
monomers to the filament tip, but treats only the diffu- 
sive behavior of the variable describing the distance be- 
tween the obstacle and the tip. It is assumed that once 
this distance exceeds a, that the monomers can enter 
with a fixed probability independent of the tip-obstacle 
distance. This assumption needs to be evaluated by ex- 
plicit treatment of the diffusion in the monomers. In 
addition, although the form of Eq. ([!]) is confirmed by 
the force-velocity relation for microtubules jll| the de- 
cay rate of the velocity with applied force was about twice 
as large as expected from Eq. ([!]). One possible expla- 
nation of this, suggested by Mogilner and Oster is 
subsidy effects between the thirteen fibers comprising a 
microtubule "protofilament" . We intend to investigate 
the extent to which other mechanisms can account for 
such discrepancies. 



II. MODEL 

Our model system contains a fiber of protein monomers 
growing perpendicular to a flat rigid obstacle in two di- 
mensions. We will be mainly interested in the actin sys- 
tem, but the basic physics of our results is relevant to 
any fiber growing against an obstacle. Our choice of two 
dimensions is dictated mainly by computational practi- 
cality: the simulations took over two weeks of CPU time 
on a Compaq 21264 processor and our preliminary stud- 
ies indicate that the three-dimensional simulations take 
about 30 times longer. The fundamental units of the 



simulation are the monomers; their internal and rota- 
tional degrees of freedom are assumed to be included in 
our effective interaction energies. The motions of the 
monomers and the obstacle are treated via Langevin dy- 
namics. The z-direction is taken as the growth axis, with 
the obstacle parallel to the x-direction. The coordinates 
of the monomer centers-of-mass are given by f^, and the 
z-coordinate of the obstacle is called Z. The Langevin 
equations for this system are n~ 1 dfi/dt = —Fi + fi(t) 
for the monomers and fi^dZo/dt = —Fo + fo(t) f° r the 
obstacle, where the /x's are mobilities, F denotes deter- 
ministic interaction forces, and fi and and fo are random 
forces satisfying 

</f (*)//(*0) = (fi'Wtf)) = 2^ 1 kT6 lJ 6(t ~ 0, (2) 
(/f(*)//(*0>=0, (3) 
and (f (t)fo(t')}=2»o 1 kT5(t-t') . (4) 

The Langevin equations are implemented with a finite 
time step At following the procedure of Ref . 1 7 1 : 



n (t + At) = n (t ) + AtmFi (t) + g(t) y/kTfu, (5) 
and Z(t + At) = Z(t) + At^ F {t) + h{t)y/kTfio , (6) 

where g(t) and h(t) are random functions with zero time 
average, satisfying 



(9 x (t)g x (t>)) = (g z (t)g*(t')) = (h(t)h(t>)) = 2AtS t 



(7) 



To implement the last set of correlations, at each time 
step we choose g x , g z , and h from a uniform random 
distribution random from — \/6Ai to \/6At. 



A. Force Laws 

The obstacle experiences an external force of magni- 
tude F axt in the — z direction. 

In the absence of reliable force fields for the monomer- 
monomer interactions, we use a simple model form for the 
interactions which has a linear filament as the lowest- 
energy structure. This form contains two-body and 
three-body interactions. The two-body interactions are 
repulsive and have the form 



V 2 (rij) = U rcp exp [-K re p(r y 



0] 



the three-body interaction energy has the form 



V^nj , f lk ) = V att exp [-KattOi 



0] 



exp [-K att (n k - a)] (a + cos%) 



(8) 



(9) 



with < a < 1. It is attractive for 9ij > cos" 1 (—a). 
The monomer-obstacle interactions have only a two-body 
repulsive component, and have the form 



W 2 { Zi ) = V ohst exp {-K ohst \ Zi - (Z - o)|) 



(10) 
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The forces are obtained as gradients of these energy 
terms. The energies are modified by subtraction of ap- 
propriate constants to force the interaction energy to go 
to zero at a cutoff distance r max (in the case of the three- 
body terms this means that the energy vanishes if either 
Tij or Tik becomes greater than 

' max ) • 

We use two parameter sets, whose values are given in 
Table Q. These two parameter sets are chosen mainly to 
sample different shapes of the "basin of attraction" for 
the addition of a monomer, and by no means exhaustively 
sample the range of possible model force fields. The first 
corresponds to a narrow basin of attraction. The large 
value of a means that the three-body terms are positive 
only for a small range of angles. This is partly compen- 
sated for by the choice of prefactors to avoid the binding 
energy becoming too small. We will call this the "hard" 
force field. The corresponding energy contours are shown 
in Figure la. The width of the basin of attraction, or the 
region over which the force pulls the next monomer into 
its minimum-energy position, is about a tenth the size of 
a monomer, which would correspond to a few A for actin 
monomers. Figure lb shows the energy contours for the 
parameters corresponding to a wider basin of attraction, 
which is about a half the size of a monomer. We call 
this the "soft" force field. For both of the force fields, 
the binding energies are very large compared to fcT, so 
that monomer subtraction from the fiber does not occur 
in the simulations. This is a reasonable approximation; 
from the measured on and off constants in Ref. 0, the 
ratio of on to off rates at physiological actin monomer 
concentrations would be less than 0.01. 

With regard to the mobilities, the only physically rel- 
evant factor is the ratio of the obstacle mobility to the 
monomer mobility, since multiplicative changes in all the 
mobilities simply serve to scale up the fiber growth ve- 
locities; these will thus factor out of our velocity results, 
which are scaled by \j\xkT. For most of our simula- 
tions, we use a mobility of the obstacle equal to that 
of the monomers for simplicity. This would correspond 
to identifying the obstacle with a part of a fluctuating 
membrane, rather than an entire rigid particle. We have 
varied the obstacle mobility in a few cases, with results 
to be discussed below. 

TABLE I. Parameters used in simulations. Energies are 
given in units of kT, ^-parameters in units of a -1 , and r max 
in units of a, the equilibrium monomer spacing. 



Force Field 


V mp 


Vatt 


Vobst 




fvatt 


^obst 


a 


''max 


Hard 


141.3 


6510 


19.14 


8.267 


4.960 


4.960 


0.940 


1.412 


Soft 


257.5 


2151 


27.44 


7.666 


4.600 


4.600 


0.770 


1.522 




FIG. 1. Energy contours for monomer approaching fiber tip 
with hard (a) and soft (b) force fields. Contour heights corre- 
spond to integer multiples of kT, with lighter corresponding to 
lower energies. The length units are nm, assuming a monomer 
step size of 2.7 nm as for actin. 



B. Filament-Growth Procedure 

A typical physiological concentration of actin (10 fjM) 
is low in the sense that the average spacing between 
actin monomers is about 60 nm, roughly 10 times the 
monomer size. This means that the probability that 
two free monomers are near enough to interact with each 
other is very small. For this reason we adopt a growth 
procedure in which only one free monomer at at time 
interacts with the tip. This is accomplished as follows. 
We start with a fiber of six monomers pointing in the z- 
direction, at their equilibrium spacing. A free monomer 
is then added at a point on a circle of radius R centered 
on the next attachment site |l8[ ] (defined as one monomer 
spacing beyond the monomer at the fiber tip) . We choose 
R = 2.5a, which places the added monomer well beyond 
the interaction range of the monomer at the tip. The rela- 
tive probabilities of monomer addition at different points 
on the circle are proportional to exp [—W2(z)/kT]. This 
weighting is accomplished by choosing a random number 
for each potential addition point; if this random number 
is less than exp [— Wi(z) jkT\, then this point is rejected 
and another one is chosen. A new point is also chosen if 
the monomer overlaps the fiber (i.e., its distance to the 
fiber is less than r max ). The system is then stepped for- 
ward in time according to the procedure described above, 
until one of two possible termination events occur: 
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• 1. The monomer diffuses outside of the i?-circle. 
In this case it is restarted on the circle as above. 
If the obstacle abuts the fiber, the position of the 
monomer is constrained to be out of the interaction 
range of the obstacle. 

• 2. The monomer attaches to the tip. In this case, 
another monomer is started on the i?-circle. 

In this way, the CPU time that is used in the simulation 
is focused on the time that the monomers spend close to 
the tip. A typical snapshot of a simulation configuration 
is shown in Figure 2. 



FIG. 
tions. 



2. Typical fiber-obstacle configuration during simula- 



velocities at other concentrations would be given by a 
linear proportionality. 
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FIG. 3. Representative plot of number of monomers in fiber 
vs. number of time steps. Obtained for Fa/kT = 1.5 and 
hard force field. 



A. Force- Velocity Relation 



One can use the computed growth rates to predict 
growth rates for low concentrations c, by simply mul- 
tiplying the computed rates by the probability P(c) of 
finding a monomer inside the i?-circle. We obtain this 
probability numerically as 

P(c) = - [ cxp(-U(^/kT)d 2 r , (11) 

c Jr<R 

where U(f) is the energy (from both fiber and obstacle) 
associated with placing a monomer at the point r, and the 
coordinates are given with respect to the next attachment 
point jlSj] . We plot our force- velocity relations in terms 
of the force acting between the obstacle and the fiber tip. 
This exceeds the external force applied on the obstacle 
by an amount corresponding to the viscous drag on the 
obstacle as it moves through the medium. The total force 
is thus given as F — F cxt + vj no- 

III. RESULTS 

Our simulations involve 10 runs, each of which involves 
the addition of 30 monomers to the fiber tip. T his cor- 
responds to a statistical uncertainty of ^/l/300 = 6% in 
the growth velocities. Typical results for the fiber length 
as a function of time are shown in Figure 3. Note that 
there are no backwards steps, because the parameters 
that are used in the force field result in an exothermic 
enthalpy for monomer addition that exceeds kT by at 
least a factor of 20. We use a concentration correspond- 
ing to one monomer per square of side 20a; in our model, 



Figure 4a shows growth velocity (solid circles) vs. ap- 
plied force, for the "hard" force field (c/. Figure la). 
For comparison, a curve proportional to exp(— Fa/kT) 
is shown. The simulation results give noticeably lower 
velocities at finite applied forces than the exponential 
prediction. The discrepancy is about 65% at Fa/kT = 
1, and 85% at Fa/kT = 2.5. The results can be 
roughly fitted to different exponential curve, of the form 
exp (—1.7 Fa/kT). Thus the growth velocity is much 
more sensitive to force than the thermal-ratchet model 
would predict. Figure 4b shows similar results for the 
soft force field (cf. Figure lb). The free-fiber growth 
velocity is about twice that for the "hard" force field, 
because the attraction basin is larger. The discrepan- 
cies between the simulation results and the analytic the- 
ory are comparable to those seen for the "hard" force 
field, but somewhat less pronounced. The discrepancy 
at Fa/kT = 2.5 is 70%, and the exponential fit curve is 
exp (—1.5 Fa/kT). The open diamonds in Fig. 4b corre- 
spond to the results of varying the mobility fio', for the 
leftmost one the mobility is doubled, and for the right- 
most one it is reduced by a factor of ten. The effects 
of these variations are very minor, as predicted by the 
"thermal-ratchet" model [ fL2| |. 
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FIG. 4. Growth rates (solid circles) for hard (a) and soft 
(b) force fields vs. total force F. Rates given in units of fikTc, 
where /i is the monomer mobility and c is the concentration. 
Force given in units of a/kT. Solid line corresponds to expo- 
nential decay (c/. Eq. (1)). Diamonds in (a) correspond to 
mobility enhanced by factor of 2 (left) and reduced by factor 
of 10 (right). Dashed curves correspond to theory of Eq. (12). 



B. Interpretation 

We believe that the discrepancies seen in Fig. 4 re- 
sult from the restriction of monomer diffusion to the 
fiber tip by the impinging obstacle. Such restriction 
will occur even when the obstacle is elevated by a dis- 
tance a or more. Figure 5 shows energy contours for a 
monomer approaching the tip, when the obstacle is ele- 
vated a distance 1.25a relative to its equilibrium position 
for Fa/kT = 1.0. The contours at at integer multiples 
of kT. The easily accessible paths corresponding to en- 
ergies less than kT are confined to a narrow band by 
the presence of the obstacle. This is expected to slow 
the diffusion to the tip. Effectively, the monomers must 
travel through a tunnel in order to get to the basin of 
attraction near the tip. Another possible explanation 
for the observed effect would be that even in the region 
with energy less than kT, there is a finite energy from 
the interaction with the obstacle. However, this energy 



is proportional to the the length scale of the interaction 
between the obstacle and the monomers. In a few cases, 
we have made this length scale five times smaller, and the 
velocities are unchanged to within a few percent. There- 
fore, this monomer-obstacle interaction energy does not 
seem to be the major factor, but rather the blocking ef- 
fects of the obstacle. 



FIG. 5. Energy contours for monomer approaching fiber 
tip with hard force field, in presence of obstacle. Contours 
are as in Fig. 1. 

To make this physical picture more precise, we have 
calculated the velocities for model fiber configurations 
in which the obstacle is held at a fixed distance from 
the fiber tip. The results are shown in Figs. 6a and b, 
for the "hard" and "soft" force fields respectively. The 
distance Z is measured in units of the monomer size, and 
the edge of the obstacle is defined as the point where 
the monomer-obstacle interaction energy is equal to kT. 
Thus when Z = 0, the interaction energy of the last 
monomer in the fiber with the obstacle is kT. In both 
cases, the velocity at Z/a = 1 is nearly zero. Only for 
Z/a > 2 is the velocity within 20% of the free-growth 
velocity. 

The appropriate generalization of Eq. (|l]) is then the 
following: 



v(F) 



v(Z)P(Z, F)dZ 



(12) 



where F is the applied force, Z is the obstacle position, 
and P(Z,F) = (const) exp (-E/kT) is the probability 
of a certain value of Z. Here the obstacle-fiber inter- 
action energy is E = Wi(z — Z) + FZ, where z is the 
z-coordinate of the last monomer in the fiber. Equa- 
tion ( |l2| ) reduces to Eq. (|l]) if v(Z) has the form of a step 
function beginning at a Z = a, and Wi is sufficiently 
short-ranged. The dashed lines in Figs. 4a and 4b cor- 
respond to a numerical evaluation of Eq. (|l2|) . For both 
force fields, the agreement with the simulation results is 
quite close, with only about 20% discrepancies occurring 
for small but non-zero forces. Thus the gradual rise of 
the velocity seen in Fig. 6, as opposed to an abrupt jump, 
is at the heart of the observed effect. 
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validity of the assumptions underlying such simulations. 
Because the bead-motion simulations include only the 
thermal energy required to achieve a certain tip-obstacle 
spacing, it is likely that the addition rate will drop off 
more rapidly with increasing force than is assumed in 
Ref. @. 

The results obtained here are also expected to have no- 
ticeable results on the structure of membranes that are 
being pushed forward by collections of actin fibers. As 
a result of random fluctuations, some fibers will eventu- 
ally get ahead of others, and these will be exerting larger 
forces on the membrane. If the velocity drops off rapidly 
with the force, then these fibers will be slowed down sig- 
nificantly. This will result in the membrane surface being 
smoother than otherwise expected. Future work should 
treat such many-fiber effects, and also explore the effects 
of fiber growth angle and branching. 
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FIG. 6. Fiber growth velocity with fixed tip-fiber spacing, 
for hard (a) and soft (b) force fields. Z is measured relative 
to point at which tip-fiber interaction energy is kT. 



IV. CONCLUSIONS 

The physics underlying the above results is general 
enough that in most systems involving fiber growth 
against an obstacle, one should expect a decay of velocity 
with applied force more rapid than the simple exponen- 
tial form (g). This may explain some of the discrepan- 
cies pointed out in connection with the measured force- 
velocity relation of Ref. pT| . However, our application 
to microtubule growth is not quantitative enough to de- 
termined whether the present effect exceeds the subsidy 
effects discussed in Ref. Q. The results obtained here 
should be useful in explaining the basic physics of motion 
based on actin polymerization. For example, in a recent 
study van Oudenaarden and Theriot [^o| have simulated 
the propulsion of plastic beads in cell extracts with a 
model based on a number of fibers exerting forces on the 
beads. In their simulations, an assumed form is taken for 
the probability of monomer addition to a fiber in terms 
of the time-averaged position of the fiber relative to the 
bead, or equivalently the force acting between the two. 
A better knowledge of the relationship between the force 
and the monomer addition rate can help pin down the 
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